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ABSTRACT 



I argue that, within the Closed-Time-Path formahsm, pinch singularities do not appear 
in truly out of equilibrium situations. 
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There are two main diagrammatic methods for the study of finite temperature field 
theory, the so caUed imaginary formaUsm (ITF) [1,2] and the real time formalism (RTF) 
[3,4,5,6,7]. The RTF has many versions (Thcrmofield Dynamics [7], Closed-Time-Path 
Formalism (CTPF) [3,4,6], Influence Functional [8]) that are equivalent to each order, 
at least perturbatively. The main advantage of the RTF over the most popular ITF is 
that it can be easily extended to deal with systems that are out of thermal equilibrium. 
The perturbative series generated by the RTF can be organised in Fcynman diagrams 
very much like the usual zero temperature field theory, the only difference being that the 
propagators acquire a 2 x 2 matrix structure having by entries Green's functions satisfying 
different boundary conditions. This immediately raises the question of the existence of 
pinch singularities: since propagators with different boundary conditions are used, what 
prevents ill defined terms like 

/I 1 
d'^P o ~~o o ~ (1) 
— + lep"^ — — le 

from appearing on a diagram ? The existence of these singularities would be a disaster 
to perturbative calculations done in the CTPF framework. Actually, an early attempt of 
studying finite temperature field theory in real time [9] (still without the matrix propaga- 
tors) failed beyond one loop because of singularities of this kind and led to the rediscovery 
of the CTPF. In equilibrium, it is known that, after the inclusion of the matrix propa- 
gators, these singularities cancel. However, it was recently pointed out [10,11] that this 
cancelation works due to the precise form of the equilibrium particle distribution function 

= • (2) 

Any other distribution function and, in particular, a time dependent one, would lead 
to uncanceled pinch singularities in the perturbative expansion. Does that mean that 
perturbation theory breaks down in non-equilibrium processes ? 

In this letter I shall argue that when physically meaningful out of equilibrium situa- 
tions are analysed, pinch singularities do not appear and demonstrate this on a example. 

To make the discussion more concrete, let us consider a model with two scalar fields 
(p and a with masses m and M that are free for times t < and at thermal equilibrium 
at temperatures T and T' respectively. At t = an interaction of the form —\(f)a"' is 
turned on. Through their mutual interaction, the fields will then thermalize. This process, 
however, takes time. At any finite time after the interaction starts, the Green's functions 
will reflect the non-equilibrium nature of the problem. For instance, since time translation 
invariance is broken, the propagator D{x, y) will depend on the time coordinates and 
separately and not only on the difference — y^- This and similar models have been 
analysed in the literature before [12,13,14]. The Lagrangian is 

L = \d^^d^^^ - ^(/.^ + \d^ad^a - - e{t)Ma^ , (3) 

where 9{t) is the step function. The Lagrangian (3) prepares the system in an initial 
condition at t = that is not in equilibirum, that is, the density matrix describing it at 
t = does not commute with the Hamiltonian for t > 0. 
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We want to find the propagator of the </> field up to order A^. In the CTPF, the free 
Lagrangian around which the perturbation expansion is built is taken as the quadratic part 
of the Lagrangian at minus infinity. Three, four,... etc., particle correlations (inexistent in 
our model) present in the initial condition as well as the interactions present at t > are 
treated as perturbations. The bare field propagator is obtained from the quadratic part 
of (3) [6,15]: 



^ ' " ^ < m<f>ix) > < T*cj>{x)<j>{y) > ) 



^ + 27rn(bo|)<^(p' - m^) 27r(^(-po) + n(bo|))<^(p' - m?) 
27r(^(po) + n{\po\mp^ - m^) - + 27rn(|po|)<^(p' - m?) ) ' 

(4) 

where (T*) T is the (anti-)time ordering operator and n(po) is the equilibrium bosonic 
distribution function at temperature T — 1/ (3 



n{po) 



1 



e/3po - 1 



(5) 



In CTPF, there are two kinds of vertices, one connecting the 1 — 1 components of the 
propagators {—iX9{t)(pa^) and another with opposite sign connecting the 2 — 2 components 
{iX9{t)(j}a'^). The connection between the bare propagator D^, the full propagator D and 
the self-energy is 



D = (D^ -1 + E) 



(6) 



We are using a notation where D{x, y) is seen as an infinite matrix. For example, 
{DD'){x,y) stands for J dzD{x, z)D' {z,y). It is convenient to change basis in (4) and 
(6) performing a similarity transformation by the matrix 



V2V-I 1 



(7) 



after which equation (6) reads 



Da 
Dr D, 



Da 
Dr D, 



+ 



Si? 





(8) 



where Da,Dr and Dc are the retarded, advanced and correlated Green's functions 



iDnix, y) =e{x^ - y^) < [(/>(x), (t>{y)] > , 
iDAix, y) =^(yO - x^) < [<l){x), (Piy)] > , 
iDc{x,y) = < [(j){x),4){y)]+ > , 



(9) 



and 
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= — — S22 

=Sii + S22 • 



(10) 



Equation (8) implies 



Dr HD'r + ^r)-' r^D%- DIT^rDI + 0{T?) , (11a) 

Da ={Dl + Sa)-' -D\- D\T.aD\ + 0(S2) , (116) 
=DrDI -^D^'^Dl -'Da- DrT.,Da 

r^Dl - DI^rDI - D%^aD\ - Dl^,D\ + 0(S2) . (11c) 



In our model, E will be given in lowest order by a graph with n u-lines joined at their 
ends (for n < 3 there are also graphs with tadpoles). Since the interaction exists only for 
positive times, S will be given by 



E(a;,y) = ^(a:°)S(a:-|/)^(yO) , (12) 

where S is the self-energy for the system in equilibrium, that means, the self-energy for 
a system with constant coupling A. The specific form of E will not be needed in our 
argument. Let us consider the perturbative expansion of Dc (He). Some insight can be 
gained looking at (He) in position space. For example, one of the terms in (11c) is 



{Dli:,D\){x,y) = J d^z d'z' D%{x - z)^(z«)Ee(z - z')D'a{z' - y)^(^'') • (13) 

Since Dr{Da) vanishes for > x^{y^ < z^')^ the integration on z'^{z'^') is bounded from 
above by x^{y^) and from below by 0. A similar thing happens with the other two terms 
in (11c), and the range of the z^, z^ integration is limited there too. We have then that Dc 
is given in first order by an integral of well behaved functions on a bounded domain: there 
are no room for divergencies (of course, the integration on z, z! may produce the usual 
ultraviolet divergencies, unrelated to the pinch singularities). Note that this would not 
be the case had the interaction existed for all times. The unlimited range of integration 
would have produced divergencies that in momentum space show up as pinch singularities. 
But fields that have interacted since t — —00 should have attained equilibrium by any 
finite time. One either has fields that started interacting at some finite time (and then 
the range of integration is limited and no singularities arise) or fields that were in contact 
since the beginning of times ( in which case they should be thermalized at any finite 
time considered and the pinch singularities cancel as shown in [10]). The situation with 
fields at two different temperatures in contact with each other for an in finite amount of 
time is clearly unphysical. The perturbative expansion of Dr and Da also always involve 
integrations on a limited domain and never contain pinch singularities. 
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These remarks will be clearer in a more explicit calculation in momentum space. In 
particular it will be clear that pinch singularities are actually related to the unboundness 
of the domain of integration. 

Up to leading order, Df. (11c) involves products of three functions of the form 



{D^i:jDl){x,y)^J d'z d'z'D^{x-z)e{z')E,{z-z')e{z'')D'i{z'-y) 

- ! dz dz' a( Ji\a( z^'\p-^P{^-z)-ik{z-z')-iq{z' -y) 

-J azaz ^27^)4 (27r)4 >^^^ 

xD^ip)i:,ik)Dliq) 

(14) 

where k = A, R, c. Equation (14) can be evaluated with the help of 



i 



dz e*"^ = lim (15) 

?7-»0+ a + IT] 



givmg 



(D^T, D?)(t v) = lim f ^ ^ ^„-ipoXo+iqoyo+ip.{x-y) 



Po-qo + 2ir] 

(16) 



The limited range of integration on z^, z^' produces the denominator in (16) instead of the 
usual 5{po — go) and the energy is not conserved in the vertices. The physical origin of 
this is clear, at the instant t = the system suddenly acquires some " potential " energy 
it did not have before the interaction was turned on. It may seem that the substitution of 
^ipo ~ Qo) by {po — Qo + 2ir])~^ is not enough to get rid of the pinch singularities since 

'—7^ = V {—^) + ^5{po - go) , (17) 

Po-qo + 2ir] \Po-qoJ 

and the term with the S function, when used in (16), produces a term proportional to 
DiipojPj^jipoiPjD^iPOj (t)j that contains pinch singularities. An explicit calculation will 
show that the term coming from the principal value in (17) is also singular in a way that 
the whole expression (16) is well defined. 

First let us calculate the term in (11c) giving by (16) with i = R^j = c and k = A 
(we are assuming from on that > the other cases can be handled similarly): 
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lim lim 



^"^P ^gO „-ipoXo+iqoyo+ip.(x—y) ^ 



Po-qo + '^irj 



X 



= lim 



74 



e^o+J (27r) 



_g-»Poa;o+i(u'p+je)j/o 



(po + ^e)2 - a;2 (go - ie)^ - 
-ipo(a;o-?/o) ! ! v ^ 

(po + ie)2 - a;2 (po - ie)2 - a;2 (I8) 



1 



Po - t^^p - «e {Po + ^e)^ - ti^p 2(a;p + ie) 
i i 1 



— Eci-ojp + ie,p) 



_|_g-»Poa;o+i€pj/o 



where ujp = a/^? + and E^(ep,|?) is the residue of Sc at its pole (in the upper half 
plane). The first term in (18) is all we would have in the case of time independent coupling. 
There is a pinch singularity there. However, the remaining terms are singular too, and the 
singularities cancel. In fact, the integration on po gives 



(2^ 



Jp.{x-y) 



4e Up —ie 

giu,^(a;o-J/o)-6(xo-2/o)l 1 tJ-Up - ie,p) 

4e u>p + ie 



4e ujp 

8a;| 
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8u>t 



1 1 



4e ujr. 



+ e 



-iojpXQ+iepyo _ 



ujp - ep - a;2 2a;p 



+ e 



+ e' 



2 ~:2~^ci^p^p) 



<^p + - 2a;p 



-i£p{xo-yo) _ 



(19) 

In (19), the first two lines correspond to the fist line in (18). Thus, in accord to the general 
argument above, the pinch singularities cancel in this expression. Note that the terms that 
depend on xq and yo separately, and not only on the difference xq — yo, come from the 
terms in (18) that would not be present if the coupling were constant. Those terms are 
the ones that break time translation invariance and correspond to truly non-equilibrium 
effects. 

The same kind of cancelation occurs in the other two terms on (11c). We just quote 
the results of the integrations here: 



iD%^nDt)ix,y) 



and 



„ip.(x-y) 



^-i{oJp-ie){xo-yo)J_f^^^ _ 



+ 



d 



dqo (go + <^p)^ 



_ ^-i^A=^o-yo)-e{xo+yo)J_f^^^^^^^ 



+ ^i^A^o-yo)-<xo+yo) + 



Aojpe' 



+ e 



iu>p{xo+yo) . 



8a;3 



{DlT.AD\){x,y) = {DlT.RDl){y,x) , 



VI riO> 



(20) 
(21) 
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where 



/(go) = (l + 2n(go))2fl(go,p). (22) 



Here too, the new terms generated by the time dependence of the coupling contain a part 
dependent only on the difference xq — yo that cancels the singularity of the " equilibrium 
" term and the result is well defined. Note that, at no point we used the specific form of 
E, so this cancelation seems to be a general feature. In particular it is independent of the 
specific distribution function of the cr-particles and its initial temperature. For n = 1 our 
model is trivially solved and the result, expanded to order agrees with (19), (20) and 
(21). 

I have argued that the perturbation series is well defined and free of pinch singularities. 
That does not mean though that perturbation theory is adequate for all non-equilibrium 
phenomena. Looking at (19) and (20) we see that Dc calculated up to order contain 
terms proportional to xq and yo- Higher order calculations of Djk^j^^^Dc would contain 
higher powers of x^,y^. This is the way perturbation theory has of approximating the 
true large time behavior of the propagator (dumped exponential) we expect on physical 
grounds. Any finite order in perturbation theory will do a poor job in approximating 
this dumped behavior at large times. This is analogous to the fact that the series 1 — 
7^ + 7^t^/2 — ... is not a good approximation to e"'''* for t >> 7~^. A solution for this 
problem is to use full propagators -Dfl(A) = {^r^a) '^R{A))~'^j with S^.^^) calculated 
in perturbatively, instead of the bare ones. This resums these powers of x^,y^ into an 
dumped exponential we expect for stable systems. One may feel tempted though, to lift the 
restriction imposed by 9{z^),9{z^ ) in , for instance, equation (13), since the contribution 
coming from z^,z^' <0 is small for large x^, y° due to the dumped exponential behavior of 
Dr(a) • This is just a consequence of the fact that for large times the system looses memory 
of the initial condition at t = (the system thermalizes) . But in doing this approximation, 
namely, integrating over z'^^z'^' from — oo to oo, one throws away all pieces that break 
time translation invariance and the result is just the thermalized (equilibrium) Green's 
functions. For short times (shorter than the thermalization time), where the essentially 
non-equilibrium behavior occurs, perturbation theory works finely, but it is important to 
keep the restriction z^,z^ > 0. This way the non-equilibrium physics (breaking of time 
translation invariance and the dependence on initial conditions) is preserved and pinch 
singularities do not appear. 

The author would like to thank A. Das for comments. 
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